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Abstract 

For zero energy, E = 0 , we derive exact, classical and quantum solutions for all power-law 
oscillators with potentials V(r) = 7 > 0 and —00 < v < 00 . When the angular mo- 

mentum is non- zero, these solutions lead to the classical orbits p{t) = [cos/i(<^>(t) — <^o( 0 )] 
with p = i//2 - 1 ^ 0. For v > 2 , the orbits are bound and go through the origin. We cal- 
culate the periods and precessions of these bound orbits, and graph a number of specific 
examples. The unbound orbits are also discussed in detail. Quantum mechanically, this sys- 
tem is also exactly solvable. We find that when v > 2 the solutions are normalizable (bound), 
as in the classical case. Further, there are normalizable discrete, yet unbound , states. They 
correspond to unbound classical particles which reach infinity in a finite time. Finally, the 
number of space dimensions of the system can determine whether or not an E — 0 state is 
bound. These and other interesting comparisons to the classical system will be discussed. 


1 Introduction 

This all really started in Moscow, in 1992. That is where I met my colleague, Jamil Daboul, at 
the Second International Workshop on Squeezed States and Uncertainty Relations. It was held 
at the Conference Center-Hotel that the Russian Academy of Sciences uses. Late at night Jamil 
and I would get into sessions on life, physics, women, politics - you know, the usual stuff - while 
we drank his scotch. 

The physics came around to musings about why certain problems can be solved exactly while 
others cannot, and the symmetries associated with such problems. There is a “folk-theorem I 
often think of, which certainly is not exact but also certainly is intriguing. This theorem declares 
that if you can solve (or not solve) something classically the same is true quantum mechanically, 
and visa versa. 
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*** 

Things stood there until Jamil visited me last year, and we took things up again. While 
wondering about the role of the Runge-Lenz vector in potential systems, a number of small 
observations started us down the line: a) the classical orbit for the attractive — 7 /r 4 potential 
with centripetal potential barrier can be solved exactly; b) this type of quantum system is usually 
only discussed for E ^ 0; c) for E = 0, both the classical and quantum equations are simpler. 
Eventually we found out that all potentials of the form: 

V ( r ) = - ’ T ,>0 - -oo<i/<oo, (1) 


can be solved exactly, both classically and quantum mechanically, for zero binding energy, E — 0. 
- In what follows, it will be useful to switch back and forth between the variables v and g related 

by 

p,-{v-2)/ 2 , — + ( 2 ) 

Therefore, contrary to the usual scenario of solving a particular potential for all energies, one 
could solve an infinite system of potentials for a particular energy. The physics that came out 
was most amusing. In this piece I will report on this work. For further details you can consult a 
letter on the new results for the quantum system [ 1 ], as well as longer articles on the classical and 
quantum physics involved [2, 3]. 

In Section 2 I will demonstrate the solution to the classical problem. Section 3 contains 
amusing specific examples of some of the classical trajectories. The quantum solution will be 
given in Section 4. Section 5 contains interesting aspects of the quantum solutions, and then I 
give a brief closing comment. 

Before continuing, I wish to further parametrize the power-law potentials as 


V(r) 



2 ma 2 p v ’ 


p = r/a . 


( 3 ) 


The dimensional coupling constant, 7 , is more useful in classical physics. The dimensionless 
coupling constant, g 2 , is more useful in quantum physics. Note, in particular, that the constant L 0 
becomes h in quantum physics. Finally, the “effective potential,” including the angular-momentum 
barrier, is 


f/(L, r) = 


y 2 

2mr 2 r y 


( 4 ) 


2 Classical Solution 


Let us now obtain the classical solution. By substituting the angular-momentum conservation 
condition 

= L/(mr 2 ) (5) 

into the energy conservation condition 


e-r-T? 



+ r 


(6) 
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one obtains 


( 7 ) 


( dr\ 2 2 2 m(E-V)r 4 

UJ +r — 

This is essentially a first-order differential equation, which could be formally integrated. How- 
ever, for E = 0 , it is much more efficient to solve Eq. (7) directly. Converting to the dimensionless 
variable p = r ja and substituting V into Eq. (7), we obtain 


p> = p^) = pV~ 2 ">. 

“¥7 


( 8 ) 


For */ = 4 the right-hand side of this equation is unity, so the solution is a cosine. This is the 
circular orbit p = cos 9 which we will discuss in the next section. Guided by this we multiply Eq. 
( 8 ) by p 2tl ~ 2 to yield 


P 


t -i <h_ 

dxp / 


+ /) 2m = 


dp? 


+ (yf 


i 


(9) 


Now, p tl satisfies the differential equation for the trigonometric functions. Therefore, the general 
solution of Eq. (9) is given by 


= cos p(p> 


' o ) 


COS 


v — 2 


[ip — 9 q ) 


(10) 


or 


p = [cos 


(i (<p - y?o)] 1/V = 


COS 


( v * 1 ' 5 ■ y 


2 

— 2 


(11) 


The phase, tp 0 , is the integration constant. 

Actually, for bound trajectories, which are the case for v > 2, the angle <p and the phase 90 
both change value at the origin. There a particle is both at the end of a particular orbit (which 
starts and ends at the origin) and also at the beginning of the next orbit. 9 changes value because 
of the use ot polar coordinates and 90 does because of the singular nature of the potentials. One 
has to be careful in matching solutions for bound orbits, and I refer you to Ref. [ 2 ] for the details. 
For now just note that this problem can be taken care of, and we set 90 = 0 for the first orbit. 


3 Classical Trajectories 


3.1 Bound trajectories: 2 < v or 1 < \i 

For 2 < v or 1 < p, the trajectories go out of and back in to the origin in a finite amount of time. 
The reason for this is that the dynamic potential dominates at the origin, but the centripetal 
barrier dominates at a finite distance. The effective potential then asymptotes to zero from above 
as r — ► co. This is shown in Figure 1. 

These bound orbits have an opening angle at the origin of 




27T 

v — 2 


7T 

P ' 


( 12 ) 
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FIG. 1. The effective potential obtained from Eq. (4) for v = 4 in units of To/2, 
as a function of p = r/a. The form is U(p ) = 4 j p 1 — 1 / p 4 . 


The precession per orbit is 




which means that if v is a rational fraction, the trajectory will close after a finite number of orbits. 
The classical period of an orbit is 


ma 2 v/T r(6) 

Tt/= [Tj Hr(i + i/ 2 ) ’ 



^ + 2 
2v -4 


> 0 . 


(14) 


(Once again, see Ref. [2] for details.) 

Starting with very large the first orbit describes a very thin petal. The second orbit precesses 
by almost —i r, being a thin petal almost on the opposite side of the first orbit. As v gets smaller, 
the petals become larger and the precession per orbit becomes smaller. 

For example, the v — 8 case, has three petals. Here a petal is 7r/3 wide and the precession per 
orbit is — 27r/3. Thus, there are three orbits before the trajectory closes. Note that here the three 
petals in a closed trajectory cover only half of the opening angle from the origin. We show this in 
Figure 2. 

The case v — 6 is very interesting. The width of a petal is 7r/2 and the precession is — 7 t/ 2 
per orbit. Here, the width of a petal and the precession are exactly such that there is no overlap 
and also no “empty angles.” It takes four orbits to close a trajectory. This is shown in Figure 3. 
We see that the physical solution consists of two perpendicular lemniscates (figure-eight curves 
composed of two opposite petals). 


16 



FIG. 2. The first three orbits for v = 8. Each orbit is precessed — 27 t/ 3 from the 
previous one, so that by the end of the 3rd orbit, the trajectory closes. In this, and 
later orbits, we show cartesian coordinates for orientation. 



FIG. 3. The first four orbits for v = 6. Each orbit is precessed — 7r/2 from the 
previous one, so that by the end of the 4th orbit, the trajectory closes. 
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When we reach v = 4, the petals have widened so much that they form a circle. The circle 
starts at the origin, travels symmetrically about the positive ;r-axis, and returns to the origin. 
The precession is zero, so the orbit continually repeats itself. In Figure 4 we show this orbit. 

As v becomes less than 4, we can think of a petal obtaining a width greater than tt, i.e., an 
orbit consists of two spirals, one out and one in, at opposite ends of the orbit. 

Consider the special case v — 3. The width of the double-spiral orbit is still given by the 
formula for 4>„, and is 2tt. Therefore, the first orbit begins and ends towards the negative .r-axis. 
The precession is 7 r, so the trajectory closes after two orbits. We show this case in Figure 5. 

As v approaches 2. the spirals become tighter and tighter and the precession (now clockwise) 
becomes larger. In fact, the spirals’ angular variation as well as the orbit’s precession both become 
infinite in magnitude as v approaches 2 . 


3.2 Unbound trajectories: v < 2 or fi < 0 

When v reaches 2, there is a singular change. First, the double spiral becomes infinite in angular 
width. But also, the joining of the two sides of the double spiral at p = l and 9 = 0 breaks down. 
It is as if a tightly-wound double spring broke. The ends spiral out to infinity. This special case 
is a Cotes’ (infinite) spiral. It takes an infinite time to reach infinity from the origin. 

When the potential parameter v just leaves that ol the infinite spiral, that is, when one barely 
has v < 2 or p < 0, there is another change. Although the two ends of the entire trajectory still 
reach to infinity and the spirals in and out almost have infinite angular widths, the distance ol 
closest approach jumps from p = 0 to p = 1 . 

As the value of u decreases, the value of the angular width ol the trajectory, now given by 
4 >„ = 7 T / 1/^ | , also decreases accordingly. By the time v = l, the angular width lias decreased to 
27 t. Eventually it becomes less than tt, meaning the orbit comes in and out in the same hall plane. 
This happens for v < 0, i.e., when the force becomes repulsive'. 

When 0 < v < 2, the repulsive centripetal barrier dominates at small r whereas the attractive 
potential V = — 7 /r*' dominates at large r. A typical shape is familiar from the Kepler problem. 
Therefore, for 0 < v < 2, the E = 0 classical orbits are all unbounded. The distance, a. now 
has a completely different interpretation. It is now the distance of closest approach. Even so, the 
formal solution (10) remains valid for negative values of p. 

As a first example consider the case v = 3/2 or p = —1/4. This orbit has a total angular 
width of 47 t. It is shown in the two drawings of Figure 6 . The large-scale first drawing shows the 
trajectory coming in from the top, performing some gyration, and going out at the bottom. The 
small-scale second drawing shows the trajectory winding around twice near the origin, with the 
distance of closest approach being one. 

A second example is the exact Kepler potential, u = 1 or p = — 1 / 2 . Eq. (10) gives 

p -1 / 2 = cos 9/2 , (15) 


so that 


= (cos 9/2 )^ 


1 + cos 9 


(16) 


This is the famous parabolic orbit for the Kepler problem with E = 0 . This orbit is shown in the 
first drawing of Figure 7. The parabola yields an angular width of 2tt, as it should. 
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If we formally set v = 0 in the expression (3), we get a negative constant potential V(r) 7 . 
Therefore, in this case the force vanishes and we have a free particle. Its orbit must be a straight 
line. However, Eq. ( 8 ) shows that one still has the same type of solution, Eq. (11). Here it is 

p = [cosy >] -1 , x = r cos (p = a . (I 7 ) 

This is the equation for a vertical straight line that crosses the ar-axis at x = a, as required by the 
initial conditions. This orbit is shown in the second drawing Figure 7, it subtending an angular 

width of 7 r from the origin. r 

For 1 / < 0 or p < — 1 the potentials V (r) in Eq. (3) are repulsive and negative- valued tor 

all r > 0, with V(r) going to -00 at large distances. Since both the potential, V(r), and the 
centripetal potential decrease monotonically, the effective potential has no minima or maxima. 
Even so, for E - 0 these unbounded orbits behave qualitatively like those for 0 < v < 2. The 
quantity a now labels the distance of closest approach and the solutions are given by the same 
expression ( 10 ), which is also valid for all fi <C 0 . 

p = [cosp<p] 1/tl = [cos , p < o . (18) 

The most famous special case of these potentials is the “inverted” harmonic-oscillator potential, 
with v = p = -2. The orbit is given by p- [cos2v ?] _1/2 , so that 

r 2 r 2 2 ■ 2 \ x2 y 2 no} 

1 = — cos2y> = — (cos>-sin <p) = — - 
a 1 d a a 

Thus, the trajectory is a special hyperbolic orbit, whose minor and major axes are equal, b = a . 
We show this orbit as the third drawing in Figure 7. Now the angular width has decreased to 7 r/ 2 . 

As the last case, we consider the orbit for v = -4 or p = -3. This orbit is shown in the last 
drawing of Figure 7. The orbit subtends an angle of tt/ 3, again as it should. One sees that as v 
becomes more and more negative, the orbits will become narrower and narrower. This is just as 
in the bound case, where the petals became narrower and narrower as v became more and more 

positive. 


4 Quantum Solution 

Consider the radial Schrodinger equation with angular- momentum quantum number /: 


ER, = 


' h 2 (<P_ 2_d_ 

2m \dr 2 r dr 


IJJ_ + 1 ) ~\ 

r 2 ) 


+ V(r) 


R, ■ 


(20) 


This Schrodinger equation is exactly solvable for the potential of Eq. (3) for all E 
oo < i/ < oo. To see this, set E = 0 in Eq. (20), change variables to p, and then 


~P 


i 2 . One finds 


0 = 


^ + '<< + *> + 4 


p‘“i 


R,(r) 


= 0 and all 
multiply by 

( 21 ) 
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This is a well-known differential equation of mathematical physics. For v ^ 2 or p / 0, the 
solution can be directly given as 



One actually has to be careful about when an absolute value of p is called for in the labels of the 
solution, and whether the J Bessel functions are called for vs. the Y function. These details are 
given in Ref. [3]. 


5 Properties of the Quantum Solution 


5.1 Normalizable bound states: 2 < v or 1 < fi 

The normalization constants for the wave functions would have to be of the form 

r°° r 2 dr T 2 f g \ 
d fill/- 1 


Nf 2 = 


j 

Jo 


p (nth) 


(23) 


Changing variables first from r to p and then from p to 2 - = gl{\p\p^ 1 ), and being careful about 
the limits of integration for all p, one obtains 


Nr 2 = -. -A 


2/n 


H Vl/d 

dz T" 

j / i+j/ 


where 

h = l 7^° ■ 

This integral is convergent and given by , 

/, = i r (g + ») 

2^ r (j + 1) r (i + ttip + 1) ’ 

if the following two conditions are satisfied: 

2/ T 1 2 

— j — j b 1 > — F 1 > 0 . 

fid 

Eqs. (26) and (27) lead to two sets of normalizable states. The first is when 


(24) 


(25) 


(26) 


(27) 


^>0 or v > 2 , / > 1/2 . (28) 

These are oidinary bound states and result because the effective potential asymptotes to zero 
from above, as in Figure 1. In this case, for E = 0, the wave function can reach infinity only by 
tunneling through an infinite forbidden region. That takes forever, and so the state is bound. 

Note that the condition on l in Eq. (28) is the minimum nonzero angular momentum allowed 
in quantum mechanics, / mjn = 1. This agrees with the classical orbit solution which is bound for 
any nonzero angular momentum. Also, the above E = 0 solutions exist for all g 2 > 0, and not 
just for discrete values of the coupling constant. 
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5.2 Free states: —2 < v < 2 or —2 < p < 0 

For — 2 < < 2 or — 2 < [i < 0 and / > 1 (as well as the solutions with l = 0 and 0 < // or 

2 < j/) the solutions are free, continuum solutions. This is in analogy to the classical case, where 
the trajectories are normal and free. 


5.3 Unbound yet normalizable states: v < -2 or // < -2 

There is another class of normalizable solutions which is quite surprising. For any / and all v < -2 
qj- n < — 2, one can verify that the conditions of Eq. (2/) are also satisfied. Thus, even though 
one here has a repulsive potential that falls off faster than the inverse-harmonic oscillator and the 
states are not bound, the solutions are normalizable ! 

The corresponding classical solutions yield infinite orbits, for which the particle needs only a 
finite time to reach infinity [2]. But it is known that a classical potential which yields trajectories 
with a finite travel time to infinity also yields a discrete spectrum in the quantum case. This 
discrete spectrum is obtained by imposing particular boundary conditions on the solutions, which 
defines a self-adjoint extension of the Hamiltonian. (See Ref. [1].) 


5.4 Bound states in arbitrary dimensions 

One can easily generalize the problem of the last section to arbitrary I) space dimensions. Doing 
so yields another surprising physical result. 

To obtain the D-dimensional analogue of Eq. (21), one simply has to replace 2 p by (D — l)p 
and /(/ + 1) by /(/ + D - 2). The solutions follow similarly as 


Ri, 


D 


i T 
i T 

^T J (^) 




W-2\pW 
9 \ 


(29) 


>1^ 


To find out which states are now normalizable one first has to change the integration measure 
from r 2 dr to r D ~ l dr and then continue as before. The end result is that if the wave functions are 
normalizable, the normalization constant is given by 






ILD 


(30) 


where 


i r | 

f \ i i 
^2 ^ tij 

1 r | 

(l+D/ 2-1 l) 

{ ImI mJ 

1 

27T 1 / 2 r 

\ ' M ) 

i r i 

(i + ^+ 

i) 


which is defined and convergent for 


21 + D- 2 


+ i > 


- + 1 


> o 


/' 


(31) 


(32) 
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This yields the surprising result that there are bound states for all v > 2 or p > 0 when 
l > 2 — D/2. Explicitly this means that the minimum allowed l for there to be zero-energy bound 
states are: 


D = 2, l min = 2, 

D = 3 , l min = 1 , 

^ 4 i lmin — 1 i 

D> 4 , l mtn = 0 . 


(33) 


This effect of dimensions is purely quantum mechanical, and exists for all central potentials. 
Classically, the number of dimensions involved in a central potential problem has no intrinsic effect 
on the dynamics. The orbit remains in two dimensions, and the problem is decided by the form of 
the effective potential, U, which contains only the angular momentum barrier and the dynamical 
potential. 

In quantum mechanics there are two places where an effect of dimension appears. The first is 
in the factor /(/ + D -+ 2) of the angular-momentum barrier. The second is more fundamental. It 
is due to the operator 


I ' qm — 


{D - 1) d 
P dp 


(34) 


This is a new contribution to the “effective potential,” and can be calculated [1]. The end result 
is that given in Eq. (33). 

The dimensional effect produces what amounts to an additional centrifugal barrier which can 
bind the wave function at the threshold, even though the expectation value of the angular mo- 
mentum vanishes. 


6 Closing Comment 

I hope you have found this discussion of anharmonic power potentials entertaining and enlight- 
ening. Jamil and I certainly have. The intuition obtained into the workings and relationships 
between classical and quantum physics has been delightful to us, to say the least. 

Thank you very much. 


References 

[1] J. Daboul and M. M. Nieto, Phys. Lett. A (submitted). 

[2] J. Daboul and M. M. Nieto, Am. J. Phys., in preparation, on the classical solutions of the 
E = 0, power-law system. 

[3] J. Daboul and M. M. Nieto, Am. J. Phys., in preparation, on the quantum solutions of the 
E = 0, power-law system. 


24 


